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Abstract 

We revisit the concavity property of the thermodynamic entropy in order to 
formulate a general proof of the minimum energy principle as well as of other 
equivalent extremum principles that are valid for thermodynamic potentials 
and corresponding Massieu functions under different constraints. The current 
derivation aims at providing a coherent formal framework for such principles 
which may be also pedagogically useful as it fully exploits and highlights the 
equivalence between different schemes. We also elucidate the consequences of 
the extremum principles for the general shape of thermodynamic potentials 
in relation to first-order phase transitions. 
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I. INTRODUCTION 



In developing the formal structure of thermodynamics, one usually starts from the max- 
imum entropy principle as the guiding principle that is used to predict the equilibrium 
conditions which apply to isolated systems. This is, actually, the way the subject was in- 
troduced by Callen in his celebrated book [1]. Crucially important in his presentation of 
equilibrium thermodynamics is the proof of the equivalence between different representa- 
tions which are based on different choices of the natural variables that can be introduced 
in order to describe the macroscopic state of the system. Such a proof implies an exten- 
sion of the extremum principle to other thermodynamic schemes. However, the approach 
that is usually pursued to justify the above equivalence, albeit physically well founded, may 
be not completely satisfactory on the formal side. More specifically, the ordinary proof of 
the minimum energy principle is formulated for a system with just one degree of freedom 
and its extension to a thermodynamic space of higher dimensionality is not, in our opinion, 
straightforward. 

This state of affairs is also probably responsible for the somewhat uncertain status of 
"thermodynamic potentials" with multiple minima that are usually introduced when dis- 
cussing first-order phase transitions [1]. Actually, it is not immediately clear how such a 
potential, which fails to fulfil the convexity requirement, is related to the fundamental equa- 
tion of the system, unless one explicitly intends to represent a Landau free energy, i.e., the 
outcome of a mean-field calculation. 

For the above reasons, we believe that it can be useful to revisit the proof of the extremum 
principles used in thermodynamics in order to place all representations on a more clear 
mathematical basis which may turn useful also for a pedagogical presentation of the subject. 

The outline of the paper is as follows: in Section 2, we provide a general proof of the 
minimum energy principle based on the concavity property of the entropy function. Then, 
upon discussing the case of a system in contact with a reservoir, we derive in Section 3 other 
forms of the extremum principle which apply to the generalized thermodynamic potentials 
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and related Massieu functions. We also analyze the pattern of singularities of, say, the Gibbs 
free energy in proximity to a first-order transition point. Some further remarks and a brief 
summary of the main results are given in the Conclusions. 

II. THE MINIMUM ENERGY PRINCIPLE 

Following Callen [1], the fundamental problem of thermodynamics is to find the equilib- 
rium state of an overall isolated macroscopic system following the removal of one or more 
internal constraints, i.e., walls restrictive with respect to the exchange of energy and, possi- 
bly, of other extensive quantities between the various parts of the system. As is well known, 
the solution to this problem can be cast in terms of the maximum entropy principle: the 
equilibrium state eventually singled out by the system is the one that maximizes the total 
entropy out of the variety of states that are compatible with the residual constraints. 

Thermodynamics essentially postulates three properties for the entropy S [1]: 1) S is 
a well-behaved, first-order homogeneous function of the extensive parameters (a property 
leading to Euler's theorem (1)); 2) S is additive over disjoint subsystems; and 3) the partial 
derivative of S with respect to the energy U is strictly positive (implying that the tem- 
perature T > 0). In particular, this latter condition allows one to express the energy as 
a function of the entropy as well as of the other extensive parameters Xi which specify 
the state of the system, in such a way that the knowledge of U (S, Xi, X 2 , . . .) is equivalent 
to that of S(U, Xi, X 2 , . . .). In the following, we shall ignore any exception to the above 
conditions such as those arising from the existence of long-ranged interactions between the 
constituent particles (additivity and, possibly, extensivity violated), or from an effective 
ergodicity breaking (which causes the unattainability of equilibrium). 

In order to set the stage for our subsequent reasoning, we consider an isolated system 
described by the energy U, the volume V, and the number of particles N as the only extensive 
parameters. For this system, the entropy is written as 

S(U, V, N) = Su(U, V, N)U + S V (U, V, N)V + S N (U, V, N)N , (1) 
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where Sx(- • •) is the partial derivative of S with respect to X and Su(U, V, N) > 0. This 
condition allows one to solve Eq. (1) in U, 

U(S, V, N) = T(S, V, N)S - P(S, V, N)V + fi(S, V, N)N , (2) 

where T(S, V, N) = U S (S,V,N) is the temperature, P(S, V, N) = -U V (S,V,N) is the 
pressure, and fj,(S, V, N) = Un(S, V, N) is the chemical potential. 

It is useful to recall that, given a function y = Y(x, a) with Y x (x, a) ^ 0, the variable x 
can be expressed, on fairly general grounds, in terms of y as x = X(y, a), with 

Y x (X(y,a),a) Y x (X(y,a),a) 
Then, the partial derivatives of S can be identified as: 

^ v > N) = ns(u t v tN)t v tN ) s (f) ^ v ' N) ; 

_ P(g(y, V, 7V), V, TV) _ / P x 

It is now possible to show that the maximum principle, along with the extensivity and 
additivity properties, underlies the concavity of the entropy function. Let A be any num- 
ber with < A < 1. Furthermore, let {U\,Vi,Ni) and (U 2 , V-i, N 2 ) identify two generic 
macroscopic states of the system. Imagine, then, to form a single isolated system by putting 
together a fraction 1 — A of (U\, V\, Aq) and a fraction A of (U 2 , V 2 , N 2 ). Once the exchange 
of U, V, and N between the two subsystems is allowed, the overall system evolves until its 
entropy reaches a value S((l — A)Lq + XU 2 , ■ ■ ■) which is larger than (or at most equal to) 
the initial value S((l — X)Ui, . . .) + S(XU2, ■ ■ •)• Thanks to the extensivity of the entropy, 
this readily implies that S is a concave function of U, V, and N: 

S((l - X)Ui + XU 2 , . . .) > (1 - X)S(U 1 , V u N l ) + XS(U 2 , V 2 , N 2 ) . (5) 

In deriving Eq. (5), we have tacitly assumed that the state space is a convex set. This 
assumption is physically reasonable as, for instance, the set U > Uo, V > 0, and A^ > is 
an open convex subset of TZ 3 [2]. 



As is well known, the occurrence of an equality sign in (5) is linked with the phenomenon 
of phase coexistence, i.e., with the occurrence of a first-order phase transition. In such a 
case, the system is macroscopically inhomogeneous and the state ((1 — X)Ui + AC/2, • • •) is 
interpreted as a mixture of the phases (U\, V±, Ni) and (C/ 2 , V 2 , N 2 ). Hence, unless two dis- 
tinct thermodynamic phases can coexist, S(U, V, N) is a strictly concave function of U, V, N 
(i.e., Eq. (5) holds as a strict inequality). 

Owing to the concavity of the entropy, the Hessian form of S is negative semidefinite (see 
Theorem 3 in Appendix A of [2]), a property that is usually expressed as d 2 S* < 0. Upon 
taking in Eq. (5) Aq = N 2 , the (U, V) Hessian of S turns out to be negative semidefinite as 
well. This condition yields the inequalities [3]: 

Suu < , S vv < , and SuuSvv - S^y > , (6) 

which are to be satisfied for all (U,V,N). These inequalities represent the conditions of 
thermodynamic stability for a system at equilibrium. It follows from Eqs. (6) that the 
const ant- volume and constant-pressure heat capacities and the isothermal and isentropic 
compressibilities are non negative quantities [1]. Similar stability conditions do also hold 
for U(S,V,N). In fact, a rather straightforward calculation along the same lines as those 
leading to Eqs. (B33)-(B35) of [2], yields: 

U SS (S, V, N) = T S (S, V, N) = -T(S, V, N) 3 Suu(U(S, V, N) , V, N) > ; 
U VV (S, V, N) = -P V (S, V, N) = -T(S, V, N) 

x [Svv(U(S, V, N), V, N) - 2S UV (U(S, V, N), V, N) ■ P(S, V, N) 
+ Suu(U(S, V, N), V, N) ■ P(S, V, iV) 2 ] > ; 
UssUvv - U 2 SV = T(S, V, iV) 4 (SuuSw - S 2 UV ) > . (7) 

As for the sign of Uyy, note that the quantity within square brackets is the value taken in 
(1, -P{S, V, AO) by the (U, V) Hessian of S relative to (U(S, V, N), V, N). 

Now, let us consider an isolated system composed of two weakly interacting subsystems 
(say, 1 and 2), not necessarily made of the same substance, which, after removing an internal 



wall, may exchange energy and one more extensive quantity (e.g., the volume V) between 
each other. Hereafter, we shall omit in the notation explicit reference to any other extensive 
parameter that is separately conserved for each subsystem. The equilibrium state eventually 
reached by the system is the state that maximizes the total entropy 

~S{U 1 , VI; U, V) = SW(U U Vl) + S®(U - U u V - VI) (8) 

with respect to the parameters of subsystem 1. In Eq. (8), U and V are the (fixed) values 
of energy and volume pertaining to the entire system. It is worth noting that the concavity 
of and implies that S as well is a concave function of Lq and V\. 

S((l- X)U[ A) + \U[ B \ (1 - A)lq (A) + XV} B) ; U, V) 
= SW((1 - \)u {A) + \u[ B \ ...) + S {2 \U - (1 - \)u[ A) - \u[ B \ . . .) 
= ^((l - \)u[ A) + \u[ B \ ...) + S {2 \(\ - \){U - u[ A) ) + \{U - u[ B) ), . . .) 

> (1 - X)S(U {A \ V} A) ; U, V) + \S{U[ B \ V} B) ; U, V) . (9) 

As a result, the (Lq, Vi) Hessian of S is negative semidefinite, which implies: 

SuiUi < , S VlVl < , and Su^Sy^ - S$ lVl > . (10) 

We now turn to the maximum condition for S. The necessary conditions for any extremal 
point (Cq , \q ) of S are: 

/ I \ (!) / 1 \ ( 2 ) 

Su^O =► (-) (U° 1 ,V 1 °) = [-j (U-UlV-V?); (11) 

S Vl = =► ^) {1 \u 1 ,V 1 ) = ^f\u-U° 1 ,V-V 1 ). (12) 

These equations merely express the well known fact that the conditions of thermal and 
mechanical equilibrium between subsystems 1 and 2 entail the same values of temperature 
and pressure for both subsystems. Furthermore, because of the concavity property, any 
extremum of S is necessarily a global maximum, which is moreover strict if S is strictly 
concave (see Theorem 1 and the corollary of Theorem 2 in Appendix A of [2] ) . 
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The solution to Eqs. (11) and (12) is generally unique (say, U?(U, V) and V?(U, V)). In 
fact, even if the final equilibrium state hosted two coexisting phases, the values of energy 
and volume of each subsystem would be uniquely determined from (11) and (12), owing to 
the fact that N\ and N 2 are fixed. Note that the derivatives of 

S(U, V) = S(U°(U, V), V?(U, V); U, V) (13) 

are well-defined only when the point of maximum of S is unique. In this case: 

Su = S Vl ■ (U°)u + S Vl ■ (V?) v + Su = Su 

= (^) ^ (U - U°(U, V),V- V?(U, V)) ee (±) (U, V)>0; 

S v = S m ■ (U°) v + S Vl ■ (V?) v + S V = S V 

= (£) ^ (U - U°(U, V),V- V?(U, V)) ee Q (U, V) . (14) 

We further notice that, if and happen to be the same function S (i.e., they pertain 
to the same substance), then S(U, V) = S(U, V, Ni + N 2 ) is a concave function of U and V . 

We now proceed to demonstrate that the maximum entropy principle can be reformu- 
lated as a minimum principle for the total energy, under a constraint on the value of the total 
entropy. To begin with, we call U^(Si,Vi) the energy function of subsystem 1, obtained 
from Si = 5W(C/i,Fi) by solving the latter with respect to U x . Similarly, let U i2) (S 2 ,V 2 ) 
be the energy of subsystem 2. The crucial step in our proof of the minimum energy prin- 
ciple will be to show that and are convex functions. To this aim, all we need to 
recall is that U^(S,V) (as well as U^) is an increasing function of its former argument, 
since ujp(S,V) = T^{S,V) > 0. Setting S A = S^(U A ,V A ) and S B = S^(U B ,V B ), the 
concavity of 

S (1 \(l - \)U A + XU B , (1 - X)V A + \V B ) > (1 - \)S (1 \U A , V A ) + \S (1 \Ub, V b ) , (15) 
can be rewritten as 

(1 - \)S A + \S B < SW((1 - \)U (1) (S A , V A ) + \UW(S B , V b ), (1 - \)V A + \V B ) . (16) 
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After evaluating U^> at the point ((1 — X)Sa + XS B: (1 — \)Va + W B ), it immediately follows 
from Eq. (16) that: 

C/«((l - X)S A + \S B , (1 - \)V A + XV B ) < (1 - A)[/«(S A , Vk) + \UW(S B , V B ) . (17) 

Under such premises, we shall now prove that, if the total entropy takes the value S = 
S(U, V), then the function 

U(S 1 , V 1 ; S, V) = (S u V,) + U^\S - S u V - V,) , (18) 

attains its minimum for V x = V?(U,V) and S x = S^(U?(U, V), V^U, V)) = S^(U,V), 
where (Ui(U,V),Vi(U,V)) is any solution to Eqs. (11) and (12). Moreover, the minimum 
value of U is U. 

We start noting, using an argument identical to that already developed for S, that U is 
a convex function of Si and V\. Therefore, in order to achieve our goal, all we need to show 
is that the first-order derivatives of U at (S®(U, V), V®(U, V)) are both zero, since then the 
convexity of U allows one to conclude that the extremum is a global minimum. 

The general expression of the first-order derivatives of U is: 

U Sl = T (1) (^, V!) - T^ 2 \S -S U V- V ± ) 

= s^(u^(s 1 ,v 1 ),v 1 )- 1 - S$(UW(S -s u v- v 1 ),v- v,)- 1 

«( ? ) (U^S-SuV-V&V-V!)-^-) (U^(S 1 ,V 1 ),V 1 ); 

(19) 

U Vl = -P {1 \S 1 , V 1 ) + pv\s -S U V- V 1 ) 

= (UW(S -SuV-V^V-Vx)- T^(S -S^V- V,) 

/ px (1) 

-(£) (U^(S 1 ,V 1 ),V 1 )-T^(S 1 ,V 1 ). (20) 
When Si = S%(U, V) and V x = V?(U, V), the energy of subsystem 1 is 

UV{S°(U, V), V?{U, V)) = U° ± (U, V) , (21) 



8 



since is the inverse of Moreover, 

S (2 \U - U°(U, V), V - V?(U, V)) = S(U, V) - S {1) (U°(U, V), V?(U, V)) 

= S-S°(U,V), (22) 

which implies: 

U i2) {S - S^(U, V),V- V?(U, V)) = U - U?(U, V) . (23) 

Given Eqs. (21) and (23), it follows from Eq. (11) that Us 1 = 0. This result, when com- 
bined with Eq. (12), yields Uy 1 = 0. Finally, the absolute minimum of U is clearly U 
(see Eqs. (21) and (23)) and this completes our proof. We further note that U(S,V) = 
mm Sl y 1 U (S 1: Vi, S, V) is the inverse function of S(U, V). In fact, for arbitrary U and V, we 
have shown that U(S(U, V), V) = U . We point out that the hypothesis according to which 
the subsystems can only exchange two extensive parameters between each other does not 
affect the generality of our proof of the minimum energy principle; rather, this restriction 
simply avoids the use of a cumbersome notation. 

A different derivation of the minimum energy principle, which does not resort to the 
convexity of U, is also viable. In this case, one must show that the (Si,V±) Hessian of U 
for V\ = V±(U,V) and Si = Si(U,V) is positive definite. Actually, what can be achieved 
this way is a weaker result, i.e., that the Hessian of U at the extremal point is positive 
semzdefinite. 

Using Eqs. (7), one immediately gets for Vi = V?(U, V) and Si = S?(t7, V): 

U VlVl = -T°(U, V) (S VlVl - 2S UlVl P°(U, V) + Su^P^U, V) 2 ) > ; 

U SlVl = -T°(U, Vf (-S UlVl P\U, V) + S UlVl ) , (24) 

where T°(U,V) = TW(S?(U,V),V?(U,V)), P°{U,V) = pW(S?(U,V),V?(U,V)), whereas 
the arguments of the second-order S derivatives are U®{U, V), V®(U, V), U, and V. Equations 
(24) also lead to: 
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Us lSl Uv lVl ~ U 2 SlVl = T°(U, Vf (SumSvm - ~S 2 UlVl ) > , (25) 

which concludes the proof that the Hessian of U is positive semidefinite in the final equilib- 
rium state. Besides the general impossibility to conclude, on account of the above inequali- 
ties, that the U extremum is a minimum (in fact, we are abstaining from using the convexity 
of U), the intrinsic limitation of the latter proof of the minimum energy principle lies in the 
fact that it only applies when the subsystems are allowed to mutually exchange at most two 
extensive parameters. In fact, only in this case the character of the Hessian of U in the final 
equilibrium state can be decided in a relatively simple way on the basis of the sign of the 
second-order derivatives. 

III. MINIMUM PRINCIPLES FOR 
OTHER THERMODYNAMIC POTENTIALS 

Thermodynamic representations other than the entropy or the energy schemes arise when 
describing the equilibrium of a system that is in contact with a reservoir. Let us consider, 
for instance, an energy reservoir (heat bath). By definition, the temperature of a heat bath 
is the same in any state, i.e., 

(S r ) Ur (U r ,V r ) = ^, (26) 

a constant number which does not depend on the energy U r or the volume V r of the reservoir. 
Hence, the entropy of a heat bath reads as 

S r (U r ,V r ) = ^ + f(V r ), (27) 

where / is an unspecified, concave function of V r . As usual, we omitted to specify the 
particle number in the notation. 

When a system with an entropy function S(U, V) is brought into contact with a heat 
bath, the joint system being isolated from the outside environment, the final equilibrium 
state maximizes the total entropy 
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S(U; U tot , V, V r ) = S(U, V) + S r (U tot - U, V r ) , (28) 

for fixed U to t = U + U r , V, and V r (we assume that a rigid and impermeable wall keeps the 
system separate from the bath). The maximum condition then reads: 

S u = ^ Su(U,V) = ±, (29) 

which is equivalent to Us(S, V) = T, U (S, V) being the inverse function of S(U, V). It might 
happen that the solution U° to Eq. (29) is not unique. However, if S(U, V) is strictly concave, 
there is a unique point of maximum U°(T, V) for S, which represents the equilibrium value 
of the system energy. In this case, the system entropy in the joint equilibrium state is also 
well-defined, being S(U°(T, V),V) = S°(T, V) and U(S°(T, V), V) — U°(T, V). 
Let us now introduce the convex function of U and V given by [4] 

F(U ;T,V) — U — TS(U, V) . (30) 

By simply looking at its derivatives, 

Fu = l-TS v and F UU = -TS UU , (31) 

it is immediately apparent that the maximum condition for S is also the minimum condition 
for F. We call F a generalized thermodynamic potential. The minimum value F(T, V) of 
F is the usual Helmholtz free energy. In fact, F(T, V) is the Legendre transform of U(S, V) 
with respect to S: 

F(T, V) = U°(T, V) - TS(U°(T, V), V) = U(S°(T, V), V) - TS°(T, V) 

= [U(S,V)-TS] s=s0(Ty) , (32) 

where we observe that S°(T, V) is the unique solution to Us(S, V) — T (a more general case 
is treated below). It is rather simple to calculate the first-order F derivatives: 

Ft = U T (T, V) - S(U°(T, V), V) - TSu(U°(T, V), V) ■ U T (T, V) 

= —S(U°(T, V), V) = —S°(T, V) ; (33) 
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F v = Uy(T, V) —T [Su(U°(T, V),V)- U V (T, V) + S V (U°(T, V), V) 

= -TS V (U°(T, V),V) = - ^of^vlv) = U v( S °( T > n V) 

= -P{S°{T, V), V) = -P°(T, V) . (34) 

In order to calculate the second-order derivatives of F we make use of Eqs. (B5) and (B7) 
of [2]: 

Frr = -^T,V) = - Uss(s J TV)v) <0 ; 

It thus follows that F is a concave function of T and a convex function of V. 

Let us now consider the case of multiple solutions to Eq. (29). For instance, it may 
happen for a particular value T c of T that Eq. (29) is solved by all U G [t^J,^], with 
C/0 = U(S%V) and U% = U(S° B ,V). This occurs if, between U° A and U%, S(U,V) is a 
linear function of U (first-order transition at temperature T c ). In this case, Us(S, V) = T c 
is satisfied for all S G [S^S^], and Us cannot be inverted as a function of S. However, 
the function F is still well-defined, along with its global minimum F(T C , V). Furthermore, 
Eq. (32) still holds, provided we call S°(T, V) the unique solution to U(S, V) = U°(T,V). 
In particular, 

lim F(T, V) = U° A - T C S° A = U° B - T C S° B = lim F(T, V) , (36) 

which means that F(T, V) is continuous for T = T c . However, F(T, V) has a cusp-like 
singularity for T = T c : 

lim F T (T, V) = -S° A ± -S° B = lim F T (T, V) . (37) 

As a further example, let us consider the case of a system exchanging energy and volume 
with a reservoir. The values of temperature and pressure are both fixed for the reservoir: 

(S r ) Ur (U r ,V r ) = ^ and (S r ) Vr (U r ,Vr) = ^ ■ (38) 
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Therefore, the bath entropy is now fully specified as 

S r (U r ,V r ) = Ur+ T PVr . (39) 
The maximum conditions for the total entropy 

S(U, V; U tot , Kot) = S(U, V) + S r (U tot - U, V tot - V) (40) 

then read: 

Su = => Su(U,V) = ^; (41) 

S y = =► Sv(U,V) = ^. (42) 

Assuming a unique solution for the above equations, one finds U = U°(T,P) and 
V = V°(T,P). Upon assuming S°(T,P) to be S(U°(T, P), V°(T, P)), one has 
U(S°(T,P),V°(T,P)) = U°(T,P). Moreover, Eqs. (41) and (42) are equivalent to 
U S (S, V) = T and U V (S, V) = -P. 

As we did previously for a system in contact with a heat bath, it is appropriate to 
introduce the auxiliary, convex function of U and V, given by 

G(U, V;T,P) = U- TS(U, V) + PV . (43) 
Clearly, G attains its minimum for U = U°(T, P) and V = V°(T, P), since 

G v = l- TS V and G v = -TS V + P (44) 

are both zero. The minimum value G(T, P) of G is the Gibbs free energy. In fact, 

G(T, P) = U°(T, P) - TS(U°(T, P), V°(T, P)) + PV°(T, P) 
= U{S°{T, P), V°(T, P)) - TS°{T, P) + PV°{T, P) 
= [U(S, V)-TS + PV} s=s0(T<nv=v0{TjP) (45) 

is the Legendre transform of U(S, V) with respect to S and V, which have been replaced by 
their conjugate variables T and — P (observe that (S°(T, P),V°(T, P)) is, by our previous 
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assumption, the unique solution to Us(S,V) = T and Uv(S,V) = —P). The first-order 
derivatives of G are simply calculated as: 

G T = C/°(T, P) - S°{T, P)-T [Su(U°(T, P), V°{T, P)) ■ C/°(T, P) 
+ S V (U°(T, P), V°(T, P)) ■ V$(T, P)] + PV$(T, P) 

= -S°(T,P); (46) 

G P = U° P (T, P)-T [Su(U°(T, P), V°(T, P)) ■ U° P (T, P) 

+ S V {U°{T, P), V°(T, P)) ■ V£(T, P)] + V°(T, P) + PVp>{T, P) 

= V°(T,P), (47) 

whereas, with the help of Eqs. (B20)-(B22) and (B24) of [2], the second-order G derivatives 
turn out to be: 

Gtt - -b T {l , P) — — — 2 < , 

UssUvv — Usv 

c - V o (TP) - Uss(S°(T,P),V°(T,P)) 
Gpp - V P (1 , f) — — — 2 < U , 

GttGpp — G 2 TP = (UssUvv — f|y) > . (48) 

Hence, G is a concave function of both T and P. 

Summing up, when a system is in thermal and mechanical contact with a reservoir, it 
is G (the generalized Gibbs potential) that is minimum at equilibrium, not the Gibbs free 
energy as is sometimes stated. Similarly, if the wall between the system and the reservoir 
is permeable to the flow of energy and particles while being restrictive to volume, it is 
A(U, N; T, V, jj) = U — TS{U,V, N) — jiN that is minimized at equilibrium, its minimum 
value being the system grand potential. 

Needless to say, equivalent maximum principles hold for the functions S(U, V) — (1/T)U 
and S(U,V) - (1/T)U - (P/T)V, equal to —(1/T)F(U;T, V) and —(1/T)G(U, V; T, P), 
respectively. Their maximum loci correspond to the usual Massieu functions. 

It is worth observing that the minimization of the generalized Gibbs potential is correct 
also in the common situation of determining the equilibrium between two systems that are 
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in contact with the same energy and volume reservoir. In fact, the maximum condition for 
the total entropy 

S = SM(U 1: V 1 ) + S^(U 2 , V 2 ) + Utot - Ul ~ U2 ~ Vl ~ ^ (49) 

can be rather obviously translated into the minimum condition for the convex function 
V 1} U 2 , V 2 ; T, P) = G^(U 1} V x ; T, P) + G {2) (U 2 , V 2 ; T, P) 

= U 1 - TSW(U U V x ) + PV X + U 2 - TSW(U 2 , V 2 ) + PV 2 , (50) 

as immediately follows from computing the partial derivatives of Eq. (50). 

In closing, we re-examine the question of the shape of the thermodynamic potentials 
for a system undergoing a discontinuous phase transition. The Gibbs free energy is the 
thermodynamic potential that is usually considered for describing the phases of matter. 
This quantity stems from G after tracing the locus of its global minimum as a function of T 
and P. When a discontinuous phase transition line is approached, a piece of ruled surface 
appears in the profile of the fundamental relation, and then also in the graph of G, which 
remains convex, even though not everywhere strictly convex. This implies a discontinuous 
evolution for the location of the absolute G minimum (not for the absolute minimum itself!), 
i.e., a jump from one valley to another as soon as the coexistence line is crossed. 

While one single minimum is the rule for G in the thermodynamic limit, this is not 
generally true for a finite system. In the framework of the statistical-mechanical foundations 
of thermodynamics, this means that the microcanonical S(U, V, N) of the finite system may 
not be everywhere concave [5]. In fact, near the would-be first-order transition point, a dip 
will usually appear in the S profile which is responsible for the phenomenon of metastability. 
In turn, G is not everywhere convex and a competition arises between two different local 
minima: while the deepest minimum characterizes the most stable phase, the other one, as 
long as it is present, will be the sign that another phase is at least metastable (this is the 
usual occurrence in mean- field treatments of first-order transitions). 

The infinite-size behavior of G is sketched in Fig. 1, where its typical profile close to a 
first-order transition point is shown. Here, G(U, V; T, P) is plotted as a function of V/N at 
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constant T, for a number of values of P across the coexistence line relative to, say, the liquid 
and the vapor of a substance (only a slice of the Gibbs surface along a locus U (V) passing 
through the actual point of minimum of G is represented in the figure). The abscissa of the 
G minimum gives the specific volume of the most stable system phase for the given T and 
P values. At the same time, the G profile will show, both as a function of P at constant T 
as well as a function of T at constant P, a cusp-like singularity of the same kind discussed 
before for the Helmholtz free energy, which is responsible for the discontinuities of volume 
and entropy, respectively, across the transition. 

IV. CONCLUSIONS 

Starting from the maximum entropy principle, we have provided a detailed derivation 
of the principle of minimum energy and of similar principles for thermodynamic potentials 
and Massieu functions, by resorting to the mathematical theory of concave many-variable 
functions. This calculation was also motivated by the fact that standard reference books on 
thermodynamics usually do not give enough information about this point. In our opinion, 
proving the interchangeability of all thermodynamic representations is a necessary prerequi- 
site that allows the interested reader to fully appreciate the elegance of the thermodynamic 
formalism. 
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Figure Captions 



1 : The figure shows how the profile of G plotted as a function of the volume V is 
modified when approaching a first-order phase transition, like that from vapor to 
liquid. As the pressure changes - at constant temperature - across the coexistence 
value -P C oex, one observes a thermodynamic-stability crossover from the vapor to the 
liquid phase (see text). Note that, with the single exception of P = P COC x, the point of 
minimum of G is regular (in particular, the tangent plane is well-defined there). 
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V/N (arb. units) 
Fig. 1 (S. Prestipino and P. V. Giaquinta, The concavity...) 
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